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IMPLICIT FUNCTION RENDERING METHOD OF NONMAN I FOL D , DIRECT 
DRAWING METHOD OF IMPLICIT FUNCTION CURVED SURFACE AND 

PROGRAMS THEREOF 

BACKGROUND OF THE INVENTION 
Technical Field of the Invention 

The present invention relates to an implicit 
function rendering method of a nonmanifold, a direct drawing 
method of an implicit function curved surface and programs 
thereof . 

Description of the Related Art 

When a surface shape of a three-dimensional object 
is rendered by a computer, parametric rendering is widely 
used, which defines the shape by a coordinate value of the 
surface of the object. As examples, free curved surface 
rendering used by a computer-aided design (CAD) , polygon 
rendering used by a virtual reality model and in an 
entertainment field, etc. can be cited. The 
paramGtric Parametric rendering has a feature that a data 
structure is intuitive and easy to understand, a shape can be 
rendered by a relatively small amount of data, and use of 
high-speed drawing hardware is easy. 

On the other hand, when a process such as 
deformation of a surface shape or a set operation is carried 
out, it is often convenient to use implicit function 
rendering^ which uses a function of an implicit form to 



indirectly render a curved surface. This rendering is used 
in fields of modeling, shape processing, physical simulation^ 
etc. This rendering form has a feature that it facilitates a 
complex process^ because of its simple data structure and 
5 suited its suitability for hierarchization or parallelization 
of the process. The par amctric Parame trie rendering and the 
implicit rendering of the a curved surface have advantages 
and disadvantages. A capability Capability of selecting the 
form in accordance with a process is preferably provided^ 
10 and, in such a case, a need arises to convert both forms. 

As prior art documents concerning the present 
invention, the following^- can be cited. 
[Non-patent document 1] 

J. Bloomenthal, Polygonization of implicit surfaces, 
15 Computer Aided Geometric Design, 5:341-355,1988. 
[Non-patent document 2] 

J. Bloomenthal, Introduction to Implicit Surface, Morgan 
Kaufmann Publishers, Inc., 1997. 
[Non-patent document 3] 
20 J. Bloomenthal and K. Ferguson, Polygonization of non- 
manifold implicit surfaces, In SIGGRAPH '95, pages 309-316, 
1995. 

[Non-patent document 4] 

M. Brady, K. Jung, H. T. Nguyen and T Nguyen, Two-phase 
25 perspective ray casting for interactive volume navigation, In 
Visualization 97, pages 243-56,1997. 
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[Non-patent document 5] 

K. Engel, M. Kraus and T. Ertl, High-quality pre-integrated 

volume rendering using hardware-accelerated pixel shading, In 

Eurographics/SIGGRAPH Workshop on Graphics Hardware '01, 
5 pages 9-16,2001. 

[Non-patent document 6] 

H. C. Hege, M. Seebas, D. Stalling and M. Zockler, A 

generalized marching cubes algorithm based on non-binary 

classifications, Technical report, Konrad-Zuse-Zentrum fur 
10 Information stechnik (ZIB) , 1997 . 

[Non-patent document 7] 

Philippe Lacroute and Marc Levoy, Fast volume rendering 

using a shear-warp factorization of the viewing 

transformation, In SIGGRAPH ? 94, pages 451-458, 1994. 
15 [Non-patent document 8] 

W. E. Lorensen and H. E. Cline, Marching cubes: a high 

resolution 3d surface reconstruction algorithm, In SIGGRAPH 

f 87, pages 163-169,1987. 

[Non-patent document 9] 
20 J. Rossignac and M. O'Connor, SGC:A dimension independent 

model for pointsets with internal structures and incomplete 

boundaries, In Geometric Modeling for Product Engineering, 

pages 145-180,1990. 

[Non-patent document 10] 
25 H. Tuy and L. Tuy, Direct 2d display of 3d objects, IEEE 

Mag, Computer Graphics and Applications, 4 (10) : 29-33, 1984 . 
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[Non-patent document 11] 

A. P. Witkin and P. S. Heckbert, Using particles to sample 
and control implicit surfaces, In SIGGRAPH f 94, pages 269- 
278,1994. 
5 [Non-patent document 12] 

Shuntaro Yamazaki, Kiwamu Kase and Katsushi Ikeuchi, 
Nonmanifold implicit surfaces based on discontinuous 
implicitization and polygonization, In Geometric Modeling and 
Processing, pages 138-146. IEEE, July 2002. 
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Many methods have been presented concerning 
conversion of implicit function rendering into parametric 
rendering. Especially, conversion into a triangle mesh can 
be carried out stably and at a high speed. However, -a— 

15 mothod m ethods for converting parametric rendering into 

implicit function rendering doco do not always provide good 
results. In the conventional implicit function rendering, a 
surface to be handled is limited to two kinds of manifolds 
because of use of a continuous real valued function for 

20 defining a function field. Thus, if there are nonmanifold 
characteristics in an entry, for example, if there arc is a 
boundary on a curved surface and a branch on a surface, there 
is a problem that these portions are lost to greatly change 
shapes . 

25 On the other hand, implicit function rendering 

itself has sufficient information to represent a curved 
surface. However, since a value to be held as data is not 
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directly connected to a geometric shape of the curved 
surface, the surface must be made visible in order to check a 
real shape . When shape designing is carried out by using the 
implicit function curved surface, however, remeshing must be 
5 carried out because of a change in a function field made by a 
process such as deformation sequentially updated during the 
process of surface mixing or deformation, and^ consequently^ 
the process takes time. Additionally, as a surface shape 
becomes more complex, the number of formed polygons is 
10 increased, creating a problem of a slow drawing speed. 

SUMMARY OF THE INVENTION 
The present invention is designed to solve the 
foregoing problems and to satisfy requirements. That is, a 

15 first object of the present invention is to provide an 

implicit function rendering method of a nonmanifold, which 
converts a nonmanifold curved surface from parametric 
rendering into implicit function rendering while correctly 
holding its shape even if there are nonmanifold 

20 characteristics in an entry, for example, there arc is a 

boundary and a branch on the curved surface, and a program 
thereof therefor . 

A second object of the present invention is to 
provide a direct drawing method of an implicit function 

25 curved surface, which can deal with updating of an implicit 
function curved surface by deformation or the like at a high 
speed, carry out high-speed drawing even if a shape is 
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complex, and thereby display an implicit function field at an 
interactive speed, and a program thereof therefor . 

According to the present invention, an implicit 
function rendering method of a nonmanifold is provided, which 
5 is characterized in that an implicit function field of a 

nonmanifold is held in a form of volume data, a value of an 
implicit function at a point between lattice points is 
decided by interpolation, and if a difference in code 
distances between two adjacent voxels to be interpolated is 
10 larger than a fixed width, no surface is formed between the 
voxels . 

By this method, when there is a boundary on a curved 
surface, it is possible to prevent formation of an 
unnecessary surface by interpolation. 
15 According to a preferred embodiment of the present 

invention, only when the following relations are all 
satisfied, 

ue(-oo), t) ... (2) 

vg [t, oo] ... (3) 
20 o<( (-u) -t) + (v-t) <aw ... (4) 

but a(>l) , 

wherein w is a space between two optional sample points; and 
u and v (u<v) are values, respectively, 
there is a surface between these two points. 
25 By this method, if a difference in code distances 

between the two adjacent voxels to be interpolated is larger 
than the fixed width, it is possible to prevent formation of 
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any surfaces between the voxels. 

Additionally, a surface position q (e[0, 1]) is 
normalized so that a value can be on a lattice point of u 
when q=0 and can be on a lattice point of v when q=l, and the 
5 position q where there is a surface is obtained by the 
following equation : 

q=(t-u) / (v-u) . (5) 

By this method, even if there is a boundary on a 
curved surface, it is possible to obtain a surface position 

10 without any mistakes. 

According to the present invention, an implicit 
function rendering method of a nonmanifold is provided, which 
is characterized in that an entered curved surface is broken 
down into curved surface patches which enable determination 

15 of a front and a back, numbers are given to the front and the 
back, respectively, to be distinguished from each other, and 
a space is classified into a plurality of regions by using 
the number of a surface of a nearest point. 

By this method, even if there is a branch on a 

20 surface, it is possible to prevent elimination of a part of 
the surface connected at a branch portion. 

According to a preferred embodiment of the present 
invention, (1) an input nonmanifold curved surface is divided 
along a branch line, broken down into curved surface patches 

25 having no branches; 

(2) numbers i are allocated to the patches in an obtained 
order, a frond and a back of each patch are distinguished 
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from each other, a number i + is given to the front, and a 
number i" is given to the back; 

(3) a space is sampled by a lattice point p— ; and 
Euclid distance d E (p) to the curved surface and 

5 number i (p) of a surface of a nearest point are allocated to 
the lattice point; 

(4) for each lattice point p, i (p n ) is investigated at six 
adjacent points p n , and groups of (i(p), i (p n ) ) where 
i(p)*i(pn) are enumerated; 

10 (5) a group of new numbers ar c is substituted for the group 

of numbers pr e pared allocated above, but if the numbers which 
are first i + and i" become the same numbers as a result of 
the substitution, no substitution is carried out for a 
combination thereof, whereby numbers are arrayed in order 

15 from 0 at the cnd after said substitution ; and 

(6) in accordance with a substitution table, a region number 
i(p) is rewritten at each lattice point p, and an implicit 
function volume function of a real value is constituted 
comprised of the obtained volume region number i (p) and the 

20 Euclid distance d E (p) to the surface at each voxel. 

By this method, the space is classified into a 
plurality of regions by using the number of the surface of 
the nearest point to enable implicit function rendering of 
the curved surface where there is a branch. 

25 A distance included in a distance i is as 

follows : 

d^efDsi, D s (i + 1)] ... (6) 
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wherein D s is a width of each divided region of a real valued 
space representing a distance; and 

in a position p of each voxel, a region distance f s (p) is 

calculated from d E (p) and i(p) by the following equation: 

5 f s (p)=min(d E , 2 B -e)+2 B i(p) ... (7), 

s(>0) is set to a minute positive real number to round down 

d E (p) so that f s (p) can be included in a half-open section of 
(6) , wherein 

d E =j(x-X) 2 +(y-Y) 2 +(z-Z) 2 

10 where the coordinate (x, y , z) is a lattice point, and the 

coordinate (X, Y, Z) is the point closest to a curved surface 
from the lattice point and wherein 2 B is a range of 
permissible region distance values in one dimension . 

By this method, even if there is a branch on the 
15 curved surface, it is possible to obtain a surface position 
without any mistakes. 

Moreover, only when the followings are all 
satisfied, 

ue(2 B i, 2 B (i + l) ... (8) 
20 ve[2 B j, 2 B (j+l)) ... (9) 

0<(u-2 B i) + (v-2 B j )<aw ... (10) 
but i, j (0<i<j<n-l) , a(>l), 

wherein w is a space between two optional sample points; 
25 and u and v (u<v) are values, respectively, there is a 
surface between these two points , and 

wherein with respect to two sample points A and B, the 
designations i, j, u, v, n, and a are defined as follows: 
i = region number of the point A, 
30 j = region number of the point B, 
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u = region distance of the point A, 
v = region distance of the point B, 

n = total number of regions in which the region code 
distance is defined, 
5 a = a parameter that makes it possible to generate a 

cured surface between the points A and B, even if the 
curved surface exists between the points A and B , and 
the points closest to the curved surface do no conform 
to each other . 

10 By this method, if a difference in region distance 

functions between the two adjacent voxels to be interpolated 
is larger than a fixed width, it is possible to prevent 
formation of any surfaces between the voxels. 

Additionally, a surface position q (e[0, 1]) is 

15 normalized so that a value can be on a lattice point of u 

when q=0 and can be on a lattice point of v when q=l, and the 
position q where there is a surface is obtained by the 
following equation : 

q-(u-2 B i) / ( (u-2 B i) + (v-2 B j ) ... (11) 

20 By this method, even if there is a boundary on a 

curved surface, it is possible to obtain a surface position 
without any mistakes. 

According to the present invention, a direct drawing 
method of an implicit function curved surface is provided, 

25 which is characterized in that a texture Tf ron t representing a 
volume value of a slice front side and a texture T back 
representing a volume value of a slide backside are used to 
interpolate and display a volume value of a region surrounded 
with the slice front side and the slice backside. 
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By this method, it is possible to reflect an 
influence of the entire volume including a region of no 
slices in a drawing result, whereby moire can be eliminated. 

According to a preferred embodiment of the present 
5 invention, intersection points between a visual line and the 
slice front side and the slice backside are calculated; and 
from a textural value t front of the slice front side and a 
textual value t bac k of the slice backside, an influence of a 
volume located on the visual line between the slices on a 
10 color and a degree of transparency observed in this position 
is calculated to be displayed on a polygon. 

By this method, it is possible to obtain a correct 
result even in drawing which uses a small number of slices. 

Additionally, a process of calculating an observed 
15 color and an observed degree of transparency from the group 
of the textural value t fr0 nt and the textural value t ba ck is 
carried out beforehand; and a result thereof is saved as a 
two-dimensional texture in a graphics card on a simplified 
chart to be referred to by using a texture combining function 
20 during drawing. 

By this method, the amount of calculation is reduced 
to enable a high-speed drawing process. 

An implicit function curved surface represented by a 
region distance field volume is converted into such a form as 
25 to be used as a 3D texture; and with respect to a group of 

optional region distances constituted of the textural values 
tfront/ tback-/ a process of calculating a color and a degree of 
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transparency observed therebetween is carried out beforehand 
to prepare a simplified chart, whereby a drawing color is 
decided. 

According to this method, a correct displaying 
result can be subjected to a high-speed drawing process. 

Furthermore, according to the present invention, a 
computer program is provided, which is characterized by 
causing a computer to execute the method of claims 1 to 3. 

By using this computer program, even updating of the 
implicit curved surface by deformation or the like can be 
dealt with at a high speed, and drawing can be carried out at 
a high speed even when a shape is complex, whereby an 
implicit function field can be displayed at an interactive 
speed . 

Other objects and advantageous features of the 
present invention will become apparent in the following 
description made with reference to the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
20 FIG. 1 is a two-dimensional schematic view showing a 

surface reconstructed by a code distance field and linear 
interpolation. 

FIGS. 2A and 2B are two-dimensional schematic views 
showing a nonmanifold rendered in implicit function by a 
25 conventional method. 

FIGS. 3A and 3B are two-dimensional schematic views 
showing a nonmanifold rendered in implicit function by a 
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method of the present invention. 

FIG. 4A is a schematic view of slice generation 
along a volume axis. 

FIG. 4B is a schematic view of slice generation 
5 vertical to a drawing screen. 

FIGS. 5A and 5B are drawing views using textures on 
two adjacent slices. 

FIG. 6A is a view of an image on a display where a 
nonmanifold model is divided into regions. 
10 FIG. 6B is a view of an image on a display of a 

corresponding pre-integration table . 

FIG. 7A is a view of an image on a display showing a 
defect when a value of a region distance field is subjected 
to linear interpolation. 
15 FIG. 7B is a view of an image on the display showing 

a case where a nearest point is taken to prevent the defect. 

FIGS. 8A, 8B, 8C are views of images on the display 
showing entered triangle mesh models (upper stage) and 
visible drawings (lower stage) of the present invention. 
20 FIG. 9 is a two-dimensional schematic view similar 

to that of FIG. 1 showing a surface reconstructed by a code 
distance field and linear interpolation. 

FIGS. 10A and 10B are two-dimensional schematic 
views similar to those of FIGS. 3A and 3B showing a 
25 nonmanifold shape rendered in implicit function by the method 
of the present invention. 

FIG. 11 is an explanatory view showing a distance 
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function and its conditions when there is a boundary on a 
surface . 

FIG. 12 is an explanatory view showing conditions 
when there is a branch on the surface, 
5 FIG. 13 is a schematic view showing a region 

distance function of the present invention. 

FIG. 14 is a schematic view showing the region 
distance function and its conditions of the present 
invention . 

10 FIG. 15A is a view of an image on a display of 

linear interpolation . 

FIG. 15B is a view of an image on a display of a 

pre-integration table corresponding to the surface where 

there is a boundary. 
15 FIG. 16A is a view of an image on the display of 

linear interpolation . 

FIG. 16B is a view of an image on a display of a 

pre-integration table corresponding to the surface where 

there are a boundary and a branch. 
20 FIGS. 17A and 17B are views of images on displays 

similar to those of FIGS. 8A to 8C showing visible drawings 

of the present invention. 



DESCRIPTION OF THE PREFERRED EMBODIMENT 
25 Hereinafter, the preferred embodiments of the 

present invention will be described with reference to the 
accompanying drawings . 
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1. The present invention provides a method which uses a 
region distance function to generate an implicit function 
volume, and uses a discontinuous function to expand a value 
thereof by interpolation, thereby rendering a parametric 
5 curved surface of a nonmanifold in an implicit function. 

There is also provided a method which applies a 
volume rendering method to an implicit function field defined 
by the method of the present invention, and directly draws a 
curved surface from the implicit function field subjected to 
10 lattice point sampling. By using this method, an implicit 

function curved surface of a manifold or a nonmanifold can be 
drawn at a high speed irrespective of complexity of a surface 
shape . 

As a specific algorithm, a method which receives a 
15 triangle mesh and outputs volume data obtained by subjecting 
an implicit function field to lattice point sampling will be 
described . 

This method can be efficiently mounted by using a 
texture combining function to program a graphics card, and a 

20 high image quality curved surface can be displayed at an 

interactive speed. The speed is constant irrespective of a 
surface shape to enable handling of manifold and nonmanifold 
curved surfaces. 

In short, main features of the present invention are 

25 the following two. 

(1) There is provided a method which converts a parametric 
curved surface into implicit function rendering. Especially, 
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an algorithm is provided which forms a volume implicit 
function field from an optional triangle mesh, and extracts 
an implicit function curved surface therefrom. 
(2) There is provided a method which draws manifold and 
5 nonmanifold implicit function curved surfaces by applying a 
volume rendering method. 
2. Related Studies 

2. 1 Implicit Function Rendering of Nonmanifold Curved 
Surface 

10 Implicit function rendering focuses on handling of 

only a manifold curved surface. However, Rossignac et al . 
have pointed out a possibility of rendering a nonmanifold 
curved surface in implicit function by using a vector-valued 
function [Non-patent document 9] . Based on this, Bloomenthal 

15 et al. have presented a method for rendering a curved surface 
as an implicit function curved surface [Non-patent document 
3] . 

However, no clear methods have been provided to 
constitute an implicit field of a form to be handled by the 

20 above method from a curved surface generally given by 

parametric rendering. According to the present invention, 
definition of a nonmanifold curved surface is made clear, and 
the curved surface included in the definition is described to 
be sufficiently practical. Then, a method is provided, which 

25 converts an optional parametric curved surface included in 
the definition into implicit function rendering. 
2. 2 Visualization of Implicit Function Curved Surface 
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The implicit function curved surface is compatible 
to visualization by ray tracing. By using this method, it is 
possible to make a surface visible by a high image quality 
[Non-patent document 2] . However, by this method, drawing at 
5 an interactive speed is difficult because a process takes 
time . 

As a method for greatly increasing a drawing speed 
by simplifying rendering, Witkin et al. have presented a 
displaying method which uses a particle system. However, 
10 since this method assumes that an implicit function enables 
differentiation on a curved surface, it cannot be applied to 
an implicit function curved surface expanded to a 
nonmanif old. 

By converting the implicit function curved surface 
15 into a triangle mesh, it is possible to achieve a high speed 
of a drawing process. Bloomenthal et al. have presented a 
method which divides an implicit function field rendering a 
curved surface into tetrahedrons, and forms each tetrahedron 
in a triangle mesh [Non-patent document 1] . This method is 
20 expanded to a certain nonmanif old [Non-patent document 3] . 

On the other hand, Lorensen et al. have presented a 
method which uses hexahedron division of an implicit function 
field to form the implicit function field in a triangle mesh 
[Non-patent document 8], and it has subsequently been 
25 expanded to a nonmanifold shape [Non-patent document 6] , 
[Non-patent document 12] . 

According to these .methods, because of a need to 
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form the mesh from the function field before the drawing 
process, it is necessary to carry out remeshing when the 
function field is changed by a process such as deformation, 
and the process takes time. Additionally, as the surface 
shape becomes more complex, the number of formed polygons is 
increased, consequently creating a problem of a slow drawing 
speed . 

According to the method of the present invention, 
since the surface is made visible, even updating of the 
implicit function field by deformation or the like can be 
dealt with at a high speed. Moreover, since the drawing can 
be carried out irrespective of complexity of the shape, the 
drawing speed is constant. 

3. Implicit Function Rendering of Nonmanifold Curved Surface 
3. 1 Implicit Function Curved Surface 

For a curved surface Sc:R 3 in a three-dimensional 
space, consideration is given to a group of a real number t 
and a real valued function f (p) which satisfies the 
following : 

P gS <=> f (p)=t • . . (1) 

However, peR 3 is a point in the three-dimensional 
space. A code f can be considered to be a function for 
returning a value in accordance with a distance to S, and S 
is an equivalent surface where a value of f is t. Definition 
by a form of a right side of the equation (1) of the curved 
surface S is referred to as "implicit function rendering", 
and S is referred to as an "implicit function curved surface" 



at this time. 

When a function field defined by the real valued 
function f is processed by a computer, f must be sampled. 
Thus, sampling is carried out at an f lattice point 
beforehand to hold a result in a form of volume data, whereby 
a process for the implicit function field can be efficiently 
carried out. According to the present invention, this volume 
data is referred to as an "implicit function volume", and the 
implicit function field is defined by a volume. If the 
implicit function field is given by a volume, a value of an 
implicit function at a point between lattice points is 
decided by interpolating a value at a sampling point. That 
is, in order to define the implicit function curved surface 
by the function value t and the implicit function volume on 
the surface, a distance function for defining an implicit 
function field and an interpolation function for deciding a 
real surface position from the sampling point are both 
necessary . 

Definition of implicit function rendering for a 
certain curve surface is not limited to one in number. 
However, in order to stably constitute the implicit function 
field even if a shape of the curved surface is complex with 
respect to a sampling width, a code distance function is 
frequently used to form the implicit function field* This 
function returns a positive code to Euclid distance to the 
curved surface when a nearest point is in a front side of the 
surface, or a negative code when it is in a backside. In the 
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implicit function volume generated by using this function, a 
value between lattice points can be obtained by linear 
interpolation from a value at a peripheral lattice point. 
FIGS . 1 and 9 show two-dimensional schematic examples of this 
5 situation. In the drawings, a function value is interpolated 
only on a side of a lattice to decide positions of surfaces, 
and such surface positions are connected to reconstruct a 
surface . 

3. 2 Nonmanifold Curved Surface 

10 In a three-dimensional space, "set S is a 2- 

manifold" means that an endlessly small neighborhood at an 
optional point on S is inphase with a two-dimensional disc. 
Hereinafter, when a term "manifold" is used, it indicates a 
2-manifold. If an implicit function field is given by linear 

15 interpolation of a sampled value, a function value can be 
differentiated at an optional point in the space. Thus, a 
neighborhood at an optional point p on the equivalent surface 
can be approximated to a very small place of a normal Vf, 
and an implicit function curved surface is a manifold at this 

20 time. 

On the other hand, "S is a nonmanifold" means that 
in a region near a certain point on the surface, only one 
coordinate axis can be taken, or there are three or more two- 
dimensional local coordinate systems to be taken. 
25 A shape handled by the present invention is a 

parametric curved surface, especially a curved surface to be 
rendered in a triangle mesh. Triangles are manifolds except 



boundaries as long as self intersection is not taken into 
consideration . Accordingly, it is only a triangle boundary'' 
where there is a point which does not satisfy a nature of the 
manifold on the curved surface, and it can be classified into 
the following two cases. 

(1) Case where there is a boundary on the surface: a case 
where a point on a certain triangle boundary is not present 
on any other triangle boundaries. FIG . 2A shows a two- 
dimensional schematic example. 

(2) Case where there is a branch on the surface: a case where 
a point on a certain triangle boundary is present on two or 
more different triangle boundaries. FIG. 2B shows a two- 
dimensional schematic example. 

However, if two different triangles in an entered 
set of triangles share a point other than a boundary, these 
triangles are divided along a straight line which includes 
the shared point to prevent intersection between the 
triangles . 

As shown in FIGS. 2A and 2B, when the conventional 
method is used to convert such nonmanifold shapes into 
implicit function curved surfaces, the shapes are greatly 
changed. Thus, in order to correctly render the nonmanifold 
in implicit function, characteristics thereof must be held in 
the implicit function field, and an algorithm must be formed 
to reconstruct it as a surface. 

3. 3 Designing of Distance Function and Interpolation 
Function 
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As described above in 3 . 1, in order to render the 
curved surface in implicit function by using the implicit 
function volume, it is necessary to define a distance 
function for constructing an implicit function field and an 
interpolation function for deciding a surface position. 
Consideration will be given to a method which properly 
defines such functions to render features of the two 
nonmanifolds described above in 3. 2 within the framework of 
an implicit function curved surface. 

First, concerning the case where there is a boundary 
on the curved surface, in the conventional implicit function 
rendering, even if a boundary portion is extended, surfaces 
which have not been present originally are formed (FIG. 2A) . 
This means that unnecessary surfaces are formed by 
interpolation even while a correct surface position is 
defined in an implicit function volume defined by using the 
conventional implicit function. Thus, by expanding the 
interpolation function to prevent formation of surfaces in 
positions where the original curved surface is not present, 
it is possible to render a surface boundary (FIG. 3A, FIG. 
10A, FIG. 11) . 

In the case of using a code distance as a distance 
function, linear interpolation has been used as an 
interpolation function to decide a surface position. Thus, 
if this is expanded and a difference in code distances 
between two adjacent voxels to be interpolated is larger than 
a fixed width, no surfaces are formed between these voxels. 
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That is, if a space between two optional sample 
points is w, and values are u, v (u<v) , only when the 
followings are all established, 
ue (-00) , t) ... (2) 
5 ve [t, 00) . . . (3) 

0<( (-u)-t)+(v-t)<aw ... (4) 
there is a surface between these two points. However, a(>l) 
is a coefficient introduced to correctly form a surface even 
when conversion of a curved surface into an implicit function 
10 field includes an error. 

If an implicit function field is given by a value 
which is obtained by subjecting a lattice point sample value 
to linear interpolation, a function value at a point in the 
space conceivably includes an error about half of a maximum 
15 lattice point space. Thus, the experiment uses a=l . 5 . 

When the above three equations are satisfied, a 
position where a surface is present is obtained by the 
following equation : 

q=(t-u)/(v-u) ... (5) 
20 However, q(e[0, 1]) is normalized so that a value of a 

surface position can be on a lattice point of u when q=0 and 
on a lattice point of v when q=l. 

Next, concerning the case where there is a branch on 
the surface, in the conventional implicit function rendering, 
25 a part of the surfaces connected at a branch potion is 

eliminated (FIG. 2B) . This is caused by an impossibility of 
correctly rendering the surface position having the branch in 
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an implicit function volume defined by using a continuous 
real valued function. 

If a curved surface is a manifold, it is possible to 
define a front and a back of the surface in the vicinity of 
5 an optional point on the curved surface. Thus, the optional 
point on the surface is present in a portion held between two 
regions, i.e., a region where a nearest point is located in 
front side and a region where it is located in a backside, 
and such region division can be rendered by using a code 

10 distance function. However, since a front and a back of the 
surface cannot be defined if there is a branch on the surface 
to form a nonmanifold, especially a peripheral region of a 
point on a branch line cannot be classified into two kinds 
from the direction of the surface of the nearest point. 

15 The entered curved surface is broken down into 

curved surface patches to enable determination of its front 
and back, and numbers are allocated to the front and the back 
to distinguish them from each other, whereby the space is 
classified into a plurality of regions by using the number of 

20 the surface of the nearest point. In this way, the curved 

surface having a branch can be rendered in implicit function 
(FIG. 3B, FIG. 10B, FIG. 12) . Hereinafter, an algorithm of 
the region division will be described. 

(1) An entered nonmanifold curved surface is divided along a 
25 branch line, and broken into curved surface patches having no 

branches . 

(2) Numbers i are allocated to the patches in an obtained 
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order. A front and a back of each patch are distinguished 
from each other, a number of the front side is set to i + , and 
a number of the backside is set to i". 

(3) A space is subjected to sampling at a lattice point p, 

5 and Euclid distance d E (p) to the curved surface and a number 
i(p) of the surface of the nearest point are allocated to the 
lattice point. 

(4) For each lattice point p, i (p n ) is investigated at six 
adjacent points p n , and a group of (i(p), i(p n )) where 

10 i(p)*i(p) is set is enumerated. 

(5) Another new number is substituted for the prepared group 
of the numbers. However, if as a result of the substitution, 
the numbers which are first i + and i" become the same 
numbers, no substitution is carried out for the combination. 

15 Lastly, numbers are arrayed in order from 0. 

(6) In accordance with a substitution table, the region 
number i (p) at each lattice point is rewritten. 

An implicit function volume of a real value is 
constituted of the volume region number i (p) and the Euclid 
.20 distance d E (p) to the surface at each voxel which are 
obtained in the aforementioned manner. Accordingly, a 
"region distance" where a code distance is expanded is 
defined. A "region" in the region distance is an expansion 
of two kinds of positive and negative "codes" in the code 
25 distance to an optional number, which is realized by dividing 
a real number space representing a distance into regions. 

When a width of a divided region is D s , a distance 
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d 1 s included in a region i is represented by the following: 
d^etDsi, D S (I + 1)) ... (6) 

Such a region distance can be calculated by the 
following process. 
5 First, D S =2 B is set to decide a size of a region from B=B max - 
[log 2 n] . However, n is the number of regions. B ma x is a 
sufficiently large number, which is given as a bit number of 
a variable representing a region distance when it is rendered 
by the computer. 

10 At this time, a region distance f s (p) is calculated 

from d E (p) and i (p) in a position p of each voxel by the 
following equation : 

f s (p)=min(d E , 2 B -e)+2 B i(p) . .. (7) 

Here, e(>0) is a minute positive real number to round 
15 down d E (p) so that f s (p) can be included in a half-open 

section of (6). In the experiment, s is set to a lattice 
point space, B max =8 . 

Since the region distance function is not a 
continuous function, its value cannot be interpolated to be 
20 linear (see FIG. 13) . Thus, if the number of regions in a 

given region distance field is n, a space between two sample 
points of the distance field is w, and values are u, v (u<v) , 
only when there is i, j (0<i<j<n-l) where the followings are 
all established, 
25 ue(2 B i, 2 B (i + l)) ... (8) 

ve[2 B j, 2 B (j+l)) ... (9) 
0<( (u-2 B i) + (v-2 B j)<otw ... (10) 
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but a(>l), there is a surface between these two points. a>l 
is a coefficient used in the equation (5). At this time, a 
position where there is a surface can be obtained by the 
following equation : 
5 q=(u-2 B i)/(u-2 B i) + (v-2 B j) ) . . . (11) 

Definition of q is similar to that of the equation (5) . FIG. 
14 is a schematic view showing the aforementioned region 
distance function and conditions thereof. 
4. Visualization of Implicit Function Curved Surface 

10 4.1 Volume Rendering Method Using Slice 

By applying a volume rendering method based on a ray 
casting method [Non-patent document 10] to an implicit 
function volume, it is possible to directly draw an implicit 
function curved surface. 

15 The ray casting is relatively costly calculation 

because of a need to sample a function value on an optional 
visual line. Lacroute et al. have presented a method which 
achieves a high drawing speed by rendering a volume in a set 
of three slices vertical to x, y, z axes (FIG. 4A) , and 

20 carries out volume sampling only on a slice plane to increase 
calculation efficiency [Non-patent document 7] . 

This method can be efficiently mounted by using a 
polygon and a texture map [Non-patent document 4]. According 
to this method, a volume slice is mapped as a transparency 2D 

25 texture on a polygon, and a set of obtained semitransparent 
polygons is displayed by a blending from the deep side in 
order with respect to a visual line. In this display 
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process, drawing acceleration of a graphics card can be used 
to carry out drawing at a high speed. Additionally, in the 
case of using the graphics card, by using a 3D texture 
function to always form polygons vertical to the visual line 
5 direction without fixing a slice direction for drawing (FIG. 
4B) , it is possible to carry out higher quality drawing in 
exchange for a slight reduction in a speed. 
4. 2 Pre-integration Method [Non-patent document 5] 

A big problem of the volume rendering using the 

10 slice is that while interpolation of a volume value is 
carried out by optional resolution on the slice, the 
interpolation is carried out only at a fixed slice space in a 
direction vertical to the slice, and thus uncoincidence in 
interpolation accuracy causes striped defects or cracks 

15 called moire to be observed. This problem can be prevented 
by changing the number of formed slices in accordance with 
necessary interpolation accuracy. However, the drawing speed 
and the number of slices are in inversely proportional 
relation to each other and, in most cases, it is not 

20 realistic to form slices where moire is negligible. 

Concerning the problem, Engel et al. have presented 
a pre-integration method which simultaneously maps not only 
the volume value of the slice but also a volume value of an 
adjacent slice on each slice, and considers an influence of a 

25 volume present between the slices to eliminate moire. 

Consideration will be given to a situation where a 
volume is rendered in a set of slices, and a slice S front 
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indicated by front is drawn in polygon. According to the 
conventional slice method, this polygon is displayed by 
mapping a texture T fron t representing a volume value on the 
slice, and thus the volume value between the slices has no 
5 influence on a drawing result. 

On the other hand, according to the pre-integration 
method, a texture T baC k on an adjacent slice indicated by back 
is mapped simultaneously with the T front/ and a volume value 
of a region surrounded with S fron t and S back is interpolated by 

10 using the two textures to be displayed. As a result, an 
influence of an entire volume including a region where no 
slice is present can be reflected in the drawing result. 
When an appearance on a certain visual line is actually 
drawn, a point of intersection between the visual line and 

15 the Sf ron t, Stack is calculated and, from a textural value t fron t 
of the front side and a textural value t back of the back side, 
an influence of a volume located on the visual line between 
the slices on a color and a degree of transparency is 
calculated to be displayed on a polygon. 

20 For example, as shown in FIG. 5B, consideration will 

be given to a situation where a surface is displayed by the 
volume rendering method using the slice. When an appearance 
on a center visual line in the drawing is considered, by the 
conventional method using the slice, the surface is not drawn 

25 because of no slices present on the surface. On the other 

hand, according to the pre-integration method, since a change 
in a voxel value between the slices can be considered, a 
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correct result can be obtained even in drawing carried out by 
using a smaller number of slices. 

From the group of the t front , the t back , a process of 
calculating an observed color and an observed degreed of 
5 transparency is carried out beforehand, and a result is saved 
in a simplified chart called a pre-integrat ion table, whereby 
the amount of calculation can be reduced. Additionally, the 
chart is saved as a two-dimensional texture in a graphics 
card, and referred to by using the texture combining function 

10 during drawing to enable a high-speed drawing process. 

In setting of the graphics card to realize the 
above, used hardware is GeForce3 GPU by nVidia Inc., and a 
Texture harder function and a Register combiners function 
which constitute the texture combining function of this 

15 graphics card are programmed by using NVparse. 

In the setting, 3D texture mapping is executed on 
the slice vertical to a screen of FIG. 4B to carry out 
drawing. The volume is irradiated with a light by fon 
shading. Accordingly, it is necessary to calculate a 

20 gradient of a scalar field given as a volume beforehand. 
When a scalar value s in each voxel and a gradient of a 
curved surface (u x , u y , u x ) are given, a 3D texture having 
pixel values (s, u y , u z , u x ) is formed therefrom to be read. 

During the drawing, a multitexture function is used 

25 to read the 3D texture in texture units 0 and 1, and to 

provide a texture coordinate so as to take out the t ba ck from 
the unit 0 and the t fron t from the unit 1. In the pre- 
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integration table, a color for a group of optional scalar 
values rendered in 2D texture and given by t fr0 ntf t baC k is 
allocated to an RGB channel, a degree of transparency is 
allocated to an A channel. 
5 However, in order to interpolate a gradient between 

the slices during irradiation calculation, a coefficient is 
allocated to the B channel to interpolate the gradient, and a 
color is actually rendered by only two channels of RG. When 
this texture is read into a unit 3, a texture coordinate for 
10 this unit is set by graphics hardware (t fron t, t bac k) , and the 
pre-integration table is referred to. 

4. 3 Visualization Method of Implicit Function Curved Surface 

In order to draw the implicit function curved 
surface rendered in the region distance field volume by the 

15 pre-integration method, the implicit function volume may be 
converted into a form to be used as a 3D texture and, for an 
optional region distance group constituted of t fr0 nt, t b a C k, a 
pre-integration table may be prepared to decide an 
interpolation color therebetween. 

20 First, concerning the implicit function volume, a 

region distance is generally given by a real number. A 
distance is quantized by an 8-bit integer (0-255) to be used 
as a texture in a current graphics card. A gradient 
necessary for shading is calculated from a gradient of an 

25 original curved surface to be similarly quantized when the 
implicit function field is formed. 

The pre-integration table is prepared as a two- 
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dimensional RGB texture of a size 256x256- Colors allocated 
to the two channels of RG can be optionally set for each 
region distance group, and a single color is given when the 
surface is uniformly drawn. In the A channel, A=0 
5 (transparent) is set when there are no surfaces between the 
slices, and A=l (opaque) is set when there is a surface. 
Lastly, as an interpolation coefficient to be allocated to 
the B channel, a quantized value of qe[0, 1] by an 8-bit 
integer which is calculated by the equation (11) is used. 

10 Consideration will be given to an example of 

rendering the nonmanifold curved surface of FIG . 6A in 
implicit function to make it visible. This model is a 
nonmanifold shape where a spherical surface is fitted into a 
plane, and a space can be classified into three regions. A 

15 corresponding pre-integration table is as shown in FIG. 6B. 
An abscissa of the table corresponds to t frO nt€[0 f 255], and 
an ordinate corresponds to t baC ke[0, 255]. Colors of points 
of the table correspond to a combination of colors of the 
surface where a portion of no colors has no surfaces and 

20 drawn by a transparent color. Color strength corresponds to 
a coefficient for gradient interpolation. FIGS. 15A, 15B and 
FIGS. 16A, 16B show examples of other pre-integration tables. 

When the implicit function curved surface of the 
nonmanifold is made visible by the pre-integration method, it 

25 is necessary to prevent linear interpolation of a region 

distance of a voxel value given as a texture. For example, 
FIG. 7A shows a result of converting the nonmanifold curved 
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surface of FIG. 16A into an implicit function volume, and 
carrying out linear interpolation of the voxel value to make 
it visible. By the linear interpolation, a distance which 
belongs to a region originally nonpresent near the region 
5 boundary is formed, and a surface which is not present 
originally appears. This defect can be solved by not 
carrying out the linear interpolation of the voxel value and, 
as a result, a correct display result can be obtained as 
shown in FIG. 7B. 

10 5. Result of Experiment 

Based on the foregoing method of the present 
invention, software for converting a triangle mesh into an 
implicit function volume, and software for displaying the 
volume in an implicit function curved surface were mounted to 

15 make experiments for a plurality of models. A computer used 
for the experiment of visualization was a PC having a CPU 
Pentium 4 1.7 GHz, a main memory of 1.0 GB, VPU GeForce 4 
Ti4 600, and VRAM of 128 MB. A volume renderer was prepared 
to carry out pre-integration by using a slice horizontal to a 

20 drawing plane as visualization software. 

Each of FIGS. 8A to 8C shows a result of converting 
a nonmanifold surface model into implicit function rendering 
and making it visible. An upper stage is a parametric- 
rendered surface model given as an entry, and a triangle mesh 

25 is used. A lower stage shows a result of converting such a 
model into an implicit function volume of a size 256 3 by 
using a region distance field, and making it visible. 
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Holes are partially bored in the models shown in 
FIGS. 8A and 8B, constituting nonmanifolds where there are 
boundaries on surfaces. According to the method of the 
present invention, the curved surface having a boundary which 
5 the conventional implicit function rendering cannot deal with 
can be rendered. However, since the implicit function field 
is sampled by a volume, it is impossible to accurately render 
characteristics of a shape equal to/lower than a voxel width, 
and especially a shape near the boundary is slightly changed. 
10 This problem can be solved by increasing a volume sampling 
density. 

FIG. 8C shows a result of an experiment for a 
surface shape constituted of a plane and a semispherical 
surface. This surface is a nonmanifold having a boundary 

15 line and a branch, but it can be understood that a boundary 
and a branch can be correctly rendered by the method of the 
present invention. A value and a gradient of a distance 
function are subjected to 8-bit quantization to carry out 
nearest point interpolation. However, since drawing colors 

20 are decided by a combination of two textural values, 

quantization noise is hardly heard in a result image as long 
as a drawing surface is not extremely large. Colors during 
drawing can be set by distinguishing the front side from the 
backside for each pair of regions where there is a surface. 

25 In this example, different colors are set for all 
combinations . 

FIGS. 17A, 17B show images on displays similar to 
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those of FIGS. 8A to 8C which are made visible according to 
the present invention. 

Table 1 shows a drawing speed when visualization is 
carried out for volume data of sizes 128 3 and 256 3 . The 
speed is changed by changing the number of formed slices, and 
many experiments show results of achieving 30 fps necessary 
for an interactive speed. 

Table 1 



Size of 
volume 


128 J 


256 J 


Number of 
slices 


32 


64 


128 


256 


32 


64 


128 


256 


FPS 


>30 


>30 


25.0 


12.3 


9.9 


5.5 


2.8 


1.5 



As described above, according to the present 
invention, by using the distance function for defining the 
implicit function curved surface and the discontinuous 
function as the interpolation function, expansion is made to 
handle the nonmanifold curved surface in implicit function 
rendering. Additionally, the method which draws such defined 
manif old/nonmanif old implicit function curved surface at a 
high speed is provided. 

The entered parametric curved surface used in the 
experiment is only a triangle mesh, but the present invention 
can be applied to a free curved surface. However, tiny 
changes are necessary in an algorithm or the like for 
dividing the free curved surface into manifold patches, and 
experiments will have to be made in the future. Furthermore, 
the implicit function rendering of the curved surface is 
effective when deformation of the shape and a set operation 
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are carried out. However, these processes can be applied to 
the rendering of the nonmanifold implicit function curved 
surface defined by the method of the present invention. 

Many limitations placed on shapes to be handled in 
5 the conventional implicit function curved surface, and big 
problems of the surface boundary and the surface branch can 
be solved by the method of the present invention. 

As in the foregoing embodiment, based on the method 
of the present invention, by mounting the software for 

10 converting the triangle mesh into the implicit function 
volume and the software for rendering the volume in the 
implicit function curved surface as the computer programs, 
even updating of the implicit function curved surface by 
deformation or the like can be dealt with at a high speed, 

15 and drawing can be carried out at a high speed even if the 

shape is complex, whereby the implicit function field can be 
displayed at an interactive speed. 

As described above, the nonmanifold implicit 
function rendering method, the direction drawing method of 

20 the implicit function curved surface and the programs of the 
present invention are highly advantageous in that if there 
are nonmanifold characteristics in an entry, for example, if 
there is a boundary or a branch on the curved surface, the 
nonmanifold curved surface can be converted from parametric 

25 rendering into implicit function rendering while the shape 

can be correctly held, even updating of the implicit function 
curved surface by deformation or the like can be dealt with 
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at a high speed, drawing can be carried out at a high speed 
even if the shape is complex, and accordingly the implicit 
function field can be displayed at an interactive speed. 

Some preferred embodiments of the present invention 
5 have been described. However, it can be understood that 
scope of rights included in the present invention is not 
limited to the embodiments. On the contrary, the scope of 
rights of the present invention includes all improvements, 
modifications and equivalents specified in appended claims. 
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